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Abstract 

In this paper, we study a Lienard system of the form x = y—F{x) , y = — x, where F[x) is 
an odd polynomial. We introduce a method that gives a sequence of algebraic approximations 
to the equation of each limit cycle of the system. This sequence seems to converge to the 
exact equation of each limit cycle. 

We obtain also a sequence of polynomials Rn{x) whose roots of odd multiplicity are related 
to the number and location of the limit cycles of the system. 
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A two-dimensional autonomous dynamical system is defined by two coupled first 
order differential equations of the form : 



x = P{x,y) , y = Q{x,y) (1) 

where P and Q are two functions of the variables x and y and the overdots denote a 
time derivative. 

Such a type of dynamical system appears very often within several branches of sci- 
ence, such as biology, chemistry, astrophysics, mechanics, electronics, fluid mechanics, 
etc B I, 

One of the most difficult problems connected with the study of system (|l|) is the 
question of the number of limit cycles. A limit cycle is an isolated closed trajectory. Iso- 
lated means that the neighboring trajectories are not closed; they spiral either toward 
or away from the limit cycle. If all neighboring trajectories aproach the limit cycle, we 
say that the limit cycle is stable or attracting. Otherwise the limit cycle is unstable or, 
in exeptional cases, half-stable. Stable limit cycles are very important in science. They 
model systems that exhibit self-sustained oscillations. In other words, these systems 
oscillate even in the absence of external periodic forcing. Of the countless examples 
that could be given, we mention only a few : the beating of a heart, chemical reactions 
that oscillate spontaneously, self-excited vibrations in bridges and airplane wings, etc. 
In each case, there is a standard oscillation of some preferred period, waveform and 
amphtude. If the system is slightly perturbated, it always returns to the standard 
cycle. Limit cycles are an inherently nonlinear phenomena; they cannot occur in linear 
systems 0, |, |, |10|, 

The first physical model to appear in the literature which can be transformed to 
a system of type (0) containing a limit cycle is due to Rayleigh The following 

equation : 

^'^ + ,flA3_^V ^ = (2) 



dt'^ \3 dt dt 



that originated in connection with a theory of the oscillation of a violin string, was 
derived by Rayleigh in 1877. 



In 1927, the dutch scientist van der Pol |14] described self-excited oscillations in 



an electrical circuit with a triode tube with resistive properties that change with the 
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current. The equation derived by van der Pol reads : 

§+e(.^-l)|+.^0 (3) 

Equations (|) and are equivalent, as can be seen by differentiating (||) with respect 
to t and putting ^ = x. 



In 1928, the french engineer A. Lienard [|l^] gave a criterion for the uniqueness of 
periodic solutions for a general class of equations, for which the van der Pol equation 
is a special case : 

Lienard tranformed to a first order system by setting ^ = z, yielding 

^ = z , | = -x-/(a:). (5) 

In fact, in his proof, Lienard used a form equivalent to (|^), obtaining through the 
change of variable z = y — F{x), where F{x) = f{T)dr : 

Equation (^) is referred to as Lienard equation and both system (|^) and @ are called 
Lienard systems. They particular case of ([|). 

In 1942, Levinson and Smith [|I^ suggested the following generalization of system 



or equivalently : 



f = l/-i^(x) , % = -9{^) (7) 



dt dt y ^ V y 



Sytems (0) and ® are equivalent to : 

§ + /w| + 3W^o (9) 

which is sometimes referred to as the generalized Lienard equation. 



In this paper, we will consider the case g{x) = x and F{x) given by an arbitrary 
odd polynomial of degree m. The fundamental problem for this type of system is the 



determination of the number of limit cycles for a given polynomial F(x) [|17| , [Tq , [19 
20, 21, 22, 24| . For m = 3, i.e. for F{x) = aix + a^x^, it has been shown in [17 
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that the system has a unique hmit cycle if aiOs < and no hmit cycle if aiOa > 0. For 
m = 5 it has been shown in that the maximum number of limit cycles is two. For 
m > 5, there are no general results about the number of limit cycles of @. 

In this paper, we present a new method that gives information about the number 
of limit cycles of @ and their location in phase space, for a given odd polynomial 
F(x). This method gives also a sequence of algebraic approximations to the cartesian 
equation of the limit cycles. 

We will explain our method through the analysis of a very well known case, the 
van der Pol equation. In this case, we have : 



F{x) = e(xV3 - x) (10) 

We propose a function h2{x, y) = + gi^2{x)y + (70,2(2;), where (71,2(2;) and (70,2(2;) are 
arbitrary functions of x. Here, the second subindex makes reference to the degree of the 
polynomial /i2 with respect to the y variable. Then we calculate h2 = {y — F(x))^ — 
This quantity is a second degree polynomial in the variable y. We will choose 
(71,2(2;) and (70,2(2;) in such a way that the coefficients of y'^ and y in h2 are zero. From 
these conditions , we obtain (71.2(2;) = ki and (70,2(2;) = + /cq, where and ki are 
arbitrary constants. As F{x) is an odd polynomial, if {x,y) is a point of the limit 
cycle of (1^), then the point (— x, —y) also belongs to this limit cycle. The equation of 
a limit cycle of (^ must be invariant by the transformation {x,y) —>■ {—x, —y). We 
want the function h2{x,y) to have this symmetry too. Thus we take ki = 0. We then 
have /12 = R2{x) = —2xF{x) = — 2ex^(x^/3 — 1). The polynomial R2{x) is even and it 
has exactly one positive root of odd multiplicity, i.e. x = y/S. 

If we integrate the function h2 along the limit cycle, we have : h2{x{t),y{t))dt = 
Jq R2{x{t))dt , where T is the period; but h2{x{t),y{t))dt = h2{x{T),y{T)) — 
h2{x{0),y{0)) = 0. Consequently, we find : R2{x{t))dt = 0. This last equality 
tells us that there cannot be any limit cycle in a region of the phase plane where R2{x) 
is of constant sign. For the van der Pol system, R2{x) has a root of odd multiplicity 
at a; = y/S , hence the maximum value of x for the limit cycle must be greater than 
V3. The curves defined by h2{x, y) = + + fco = are closed for k^ < 0. 

As the next step of our procedure, we propose a fourth degree polynomial in y 
for the function hi{x,y), i.e. hi{x,y) = y^ + g3,4.{x)y^ + g2A{x)y'^ + gi,4{x)y + 5-0,4(2;) 
(polynomials hn{x,y) with n odd do not give useful information about the limit cycles 
of the system since the level curves hn{x,y) are open and the polynomials Rn{x) have 
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always a single root of odd multiplicity at x = 0) . By imposing the condition that 
must be a function of only x, we find = R^lx), where R^lx) is an even polynomial 
of tenth degree. The roots of i?4(x) depend of e, hence in the following, we will take 
e = 1. For this case, R^lx) has only one positive root of odd multiplicity, given by 
X ~ 1.824. This root is greater than the root of R2{x). Obviously, the maximum value 
of X for the limit cycle must be greater than this value. 

We have in this way a new lower bound for the maximum value of x on the limit 
cycle. Moreover the number of positive roots of odd multiplicity is equal to the number 
of limit cycles of the system. The condition that must be a function only of x imposes 
a first order trivial differential equation for each function gj n{x). These equations can 
be solved by direct integration and we obtain in this way all the functions gj^n{x). We 
take all the integration constants, that appear when we solve these equations, equal to 
zero. In this way, the level curves hi{x,y) = K are all closed for positive values of K 
and even values of n. Moreover, the function h^lx, y) is a polynomial in x and y. 

We have found the same results for greater values of n even. We have calculated 
hn{x,y) and Rn{x) up to order 20. In all cases, the polynomials Rn{x) have only 
one positive root of odd multiplicity. Let r„ be the number of such roots. For the 
van der Pol equation, it seems that r„ = 1 Vn even. These roots approach in a 
monotonous fashion the maximum value of x on the limit cycle. The functions hn{x, y) 
are polynomials in x and y for all n. The level curves hn{x,y) = K are all closed 
for positive values of K. By imposing the condition that the maximum value of x on 
the curve hn{x,y) = K > must be equal to the root of Rn{x), we find a particular 
value of K for each n even. Let us call this value K*. The level curve hn{x,y) = K* 
represents an algebraic approximation to the limit cycle. 

In fig. |I] and ^ we show this curve for the values n = 6 and n = 18, respectively. In 
table ^ we give the values of the roots of Rn{x) and the values of K* for 2 < n < 20. 
The numerical value of the maximum of x on the limit cycle, determined from a nu- 
merical integration of (P), with F{x) defined by (|ToD , is Xmax — 2.01 (e = 1). It is clear 
that the roots of Rn{x) seem to converge to Xmax and the curves /„(x, y) = K* seem 
to converge to the limit cycle. 



We have also studied the case : 
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F{x) = 0.8x- ^x^ + 0.32x^ (11) 
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This system has exactly two hmit cycles [|TB|. We have calculated the polynomials 
hn{x,y) and Rn{x) up to n = 16. The polynomials Rn{x) have exactly two positive 
roots of odd multiplicity. We conjecture that r„ = 2 Vn even. For each value of 
n, we determine two values K*-^ and K*2- The closed curves hn{x,y) = K*^ and 
hn{x,y) = K*2 provide algebraic approximations to each cycle for each value of n 
even. 

In fig. ^ and ^ we show these curves for n = 6 and n = 14, respectively. We also 
show the limit cycles obtained by numerical integration. In table 0, we give the values 
of the roots of Rn{x) and the values of K*^ and K*2 for 2 < n < 16. These roots 
seem to converge to the maximum values of x for each cycle (the numerical values 
of the maximum of x on each limit cycle are Xmax,i — 1.0034 and Xmax,2 — 1.9992 
respectively). The curves hn{x,y) = K*^ and hn{x,y) = K*2 seem to converge to each 
one of the limit cycles of the system. 

For all the cases that we have studied, we have found that the values of the constants 
K* go to zero or infinity when n ^ oo. In fact, it is easy to see from table and table 
H that the asymptotic behaviour of K* with n (for a given limit cycle), is given by 

K - a{x^axT (12) 
where a is a constant which depends on the cycle (see fig. 

We have also considered system (^ with : 

F{x) = x^ - I2x^ + X (13) 



where is an arbitrary parameter. It has been proved in |^ that this system has 
exactly two limit cycles for > 2.5. It is clear that this system has no limit cycle 
for fi < 2 because r2 = in that case. Hence, between fi = 2 and fi = 2.5 there is 
a bifurcation value fi* such that for fi < fi* the system has no limit cycles and for 
fi > fi* the system has exactly two limit cycles. When fi = fi* the system undergoes a 
saddle-node bifurcation. 

By applying our method, we can obtain lower bounds for the value of fi*. For each 
even value of n we calculate the maximum value of fi for which r„ is zero. This value 
of fi represents a lower bound for fi*. The results of these calculations are given in 
table ^ The values of /x* seem to converge very quickly, in a monotonous way, when 
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n —* oo. Numerical integrations of system with F{x) given by (|TB|) seem to confirm 

that hm„^oo fJ'n = Z^*- 

Let us point out that it is the first time, in our knowledge, that a bifurcation value 
of this type can be estimated in such a way, that is by employing an analytical method 
instead of a numerical integration of the system. 

We have also analysed system (|]) with F{x) given by : 

F{x) =x{x^ -l.6'^){x^ ~A){x^ -9) (14) 

For this case we have r2 = = 3. However, the second positive root of R^lx) is 
smaller than the second positive root of R2{x). Indeed for n = 6 we find rg = 1. An 
annihilation of two roots has occured and this phenomenon has been annonced by the 
lowering of the value of one of the roots of Rn{x). We conjecture that r„ = 1 Vra even, 
greater than 4. The numerical analysis of this system seems to indicate that it has 
exactly one limit cycle. 

For all the cases that we have studied, we have found that two types of behaviour 
of Tn are possible : 

i Tn = r'^ for arbitrary even values of n and n'. In this case the number of limit 
cycles of the system is given by this common value of the number of positive 
roots of odd multiphcity of -R„(x). 

ii the values of r„ changes with n; in this case the values of r„ decreases with n; 
moreover we have Vn — r'^ = 2p for n' > n and p G N. The roots of Rn{x) seem 
to disappear by pairs, when n increases. 

Guided by the particular cases that we have analysed, we establish the following con- 
jecture : 

Conjecture : Let be I the number of limit cycles of Let be r„ the number of 
positive roots of Rn{x) (with n even) of odd multiplicity. Then we have : 

i I < Tn^n even 

ii if n' > n then Tn — r'^ = 2p with p eN. 

We have also analysed the roots of the polynomials gj^n, with < j < n — 1. For 
odd values of j, the roots of these polynomials are also related to the number and 
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location of the limit cycles of the system. For instance, for the van der Pol equation, 
the polynomials gj,n{x) with j odd have exactly one positive root of odd multiplicity. 
These roots are an upper bound to Xmax- For a given odd value of j, the sequence of 
roots of (7j,n(x) decreases monotonously with n and seems to converge to the value of 
Xmax- The best upper bounds are given by the roots of as can be seen in table 

^. The reasons of such a behaviour of the roots of the polynomials gj,n{x) with j odd 
are not clear to us. 

We have shown in this paper that the polynomials hn{x, y) = + gn~i,n{x)y"'^^ + 
gn-2,n{x)y^~'^ + ... + gi^n{x)y + go,n{x) give a lot of information about the number and 
location of the limit cycles of (|^), in the case where F{x) is an odd polynomial (for 
the case where F{x) is not an odd polynomial, the limit cycles are not invariant under 
the transformation {x,y) {—x, —y) and the results are not conclusive). The curves 
hn{x,y) = K* give algebraic approximations to each limit cycle . These algebraic 
approximations seem to converge to the limit cycles of the system. The positive roots 
of odd multiplicity of the polynomials Rn{x) = hn{x,y) are related to the number of 
limit cycles of @ and they give lower bounds for the values of Xmax of each limit cycle. 
Moreover, the roots of gj^n{x), with odd values of j, are also related to the number of 
limit cycles and they give upper bounds to the value of Xmax for each limit cycle. 

All the relevant information about the limit cycles of seems to be contained 
in the polynomials hn{x,y). These polynomials are very easy to calculate with an 
algebraic manipulator program. 
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n 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 


root 


1.732 


1.824 


1.869 


1.896 


1.914 


1.927 


1.937 


1.944 


1.950 


1.955 


K* 

n 


3 


12.3 


54.5 


247.6 


1141 


5305 


24773 


116050 


544800 


~ 2 ■ 10*^ 



Table 1: For each value of n we give the value of the root of Rn{x) and the value of 
K* for the van der Pol equation. 



n 


root one 




root two 




2 


0.852 


0.726 


1.854 


3.439 


4 


0.905 


0.711 


1.885 


14.5 


6 


0.931 


0.739 


1.905 


67.59 


8 


0.945 


0.784 


1.920 


334 


10 


0.955 


0.840 


1.931 


1712 


12 


0.962 


0.903 


1.938 


8973 


14 


0.967 


0.974 


1.945 


47741 


16 


0.971 


1.052 


1.950 


254400 



Table 2: For each value of n, we give the two roots of Rn{x) and the values of K*^ and 
K*2 for equations (||), with F{x) given by (ITll) 
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n 


2 


4 


6 


8 


10 


12 


14 


16 


18 


20 




2 


2.057 


2.079 


2.090 


2.096 


2.100 


2.103 


2.105 


2.106 


2.107 



Table 3: We give in this table, for each even value of n between 2 and 20, a lower 
bound /i* for the value of /i*. This sequence seems to converge rapidly toward /i*. 



n 


Root of Rn 


Root of gi^n 


Root of gs^n 


2 


1.7321 






4 


1.8248 


2.2361 




6 


1.8697 


2.1924 


2.2361 


8 


1.8965 


2.1658 


2.2063 


10 


1.9144 


2.1475 


2.1854 


12 


1.9273 


2.1341 


2.1697 


14 


1.937 


2.1236 


2.1574 


16 


1.9446 


2.1152 


2.1474 


18 


1.9507 


2.1083 


2.1391 


20 


1.9558 


2.1025 


2.1321 



Table 4: For each even value of n, between 2 and 20, we give the roots of the polynomials 
Rn, gi,n and g^^n respectively, for the van der Pol equation. 
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Figure 1: The limit cycle of the van der Pol equation (exterior curve) and the algebraic 
approximation hQ{x,y) = Kq. 




Figure 2: The limit cycle of the van der Pol equation (exterior curve) and the algebraic 
approximation his{x,y) = K^^ . 
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Figure 3: The limit cycles of equation with F{x) given by ( ]Tl| ) (rough curves) and 
their algebraic approximations (smooth curves) : hQ{x,y) = K^^ and h^^x.y) = 




Figure 4: The limit cycles of equation (^ with F{x) given by ([TTl) (rough curves) and 
their algebraic approximations (smooth curves) : hi^^x, y) = K^4i and hi^^x, y) = -ft'142 
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Log[Kn] 



5 10 15 20 25 



Figure 5: We show, for the van der Pol equation, the curve Log{K*) in function of n 
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